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^ ■ ABSTRACT 

o : 

O ; 

CN I We combine COBE-DMR measurements of cosmic microwave background (CMB) 

Ch ' anisotropy with a recent measurement of the mass power spectrum at redshift z = 2.5 

1—5 ' from Lya forest data to derive constraints on cosmological parameters and test the 

• inflationary cold dark matter (CDM) scenario of structure formation. By treating 

■ the inflationary spectral index n as a free parameter, we are able to find successful 
^SJ ■ fits to the COBE and Lya forest constraints in 0,^ = 1 models with and without 

, massive neutrinos and in low-O^ models with and without a cosmological constant. 

00 I Within each class of model, the combination of COBE and the Lya forest P{k) 

I constrains a parameter combination of the form ^Imh^n^^l, with different indices 

' for each case. This new constraint breaks some of the degeneracies in cosmological 

■ parameter determinations from other measurements of large scale structure and CMB 

■ anisotropy. The Lya forest P{k) provides the first measurement of the slope of the 
Q^. linear mass power spectrum on ~ Mpc scales, = —2.25 it 0.18, and it confirms a 

basic prediction of the infiationary CDM scenario: an approximately scale-invariant 
spectrum of primeval fluctuations (n ~ 1) modulated by a transfer function that bends 
^ I P{k) towards k^~^ on small scales. Considering additional observational data, we 

^ ' find that COBE-normalized, = 1 models that match the Lya forest P{k) do not 

■ match the observed masses of rich galaxy clusters and that low-Jim models with a 
j_j ■ cosmological constant provide the best overall fit to the available data, even without 

. . the direct evidence for cosmic acceleration from Type la supernovae. With our fiducial 

parameter choices, the fiat, low-O^ models that match COBE and the Lya forest P{k) 
also match recent measurements of small scale CMB anisotropy. Modest improvements 
in the Lya forest P{k) measurement could greatly restrict the allowable region of 
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parameter space for CDM models, constrain the contribution of tensor fluctuations to 
CMB anisotropy, and achieve a more stringent test of the current consensus model of 
structure formation. 



1. Introduction 



Cosmological models based on cold dark matter (CDM) and simple versions of inflation 
have had considerable success in accounting for the origin of cosmic structure. In this class of 
models, the primordial density fluctuations are Gaussian distributed, and the shape of their power 
spectrum is determined by a small number of physical parameters that describe the inflationary 
fluctuations themselves and the material contents of the universe. For specified cosmological 
parameters, the measurement of cosmic microwave background (CMB) anisotropies by the 
COBE-DMR experiment ( Smoot et al. 199^ ; Bennett et al. 1996| ) fixes the amplitude of the 



matter power spectrum on large scales with an uncertainty of ~ 20% (e.g., Bunn &: White 1997] ) 



In this paper, we combine the COBE normalization with a recent measurement of the matter 
power spectrum by Croft et al. (1999b, hereafter CWPHK) to test the inflation+CDM scenario 
and constrain its physical parameters. A modified version of the method developed here is applied 
to a more recent power spectrum measurement by Croft et al. (2001). 

CWPHK infer the mass power spectrum P{k) from measurements of Lya forest absorption in 
the light of background quasars, at a mean absorption redshift z ~ 2.5. The method, introduced 
by Croft et al. (1998), is based on the physical picture of the Lya forest that has emerged in recent 
years from 3-dimensional, hydrodynamic cosmological simulations and related analytic models 



(e.g., Cen et al. 1994 ; Zhang, Anninos, fc Norman 1995 ; Hernquist et al. 1996 ; Bi &: Davidsen 



1997| ; plui, Gnedin, fc Zhang 19971) . By focusing on the absorption from diffuse intergalactic gas in 
mildly non-linear structures, this method sidesteps the complicated theoretical problem of biased 
galaxy formation; it directly estimates the linear theory mass power spectrum (over a limited range 
of scales) under the assumption of Gaussian initial conditions. Because the observational units 
are km s~^, the CWPHK measurement probes somewhat different comoving scales for different 
cosmological parameters: A = 2'iT/k = 2 — 12/i^^Mpc for ^lm=l, A = 3 — 16/i~^Mpc for $7^=0.55 
and f^A=0, and A = 4 - 22/i-iMpc for 17^=0.4 and S7a=0.6 {h = Hq/WO km s"^ Mpc"^). 
CWPHK determine the logarithmic slope of P{k) on these scales with an uncertainty ~ 0.2 
and the amplitude with an uncertainty ~ 35%. The extensive tests on simulations in Croft et 
al. (1998) and CWPHK suggest that the statistical uncertainties quoted here dominate over 
systematic errors in the method itself, though the measurement does depend on the assumption of 
Gaussian primordial fluctuations and on the broad physical picture of the Lya forest described in 
the references above. For brevity, we will usually refer to the CWPHK determination of the mass 
power spectrum as "the Lya P{k)." 



In the next Section, we discuss our choice of the parameter space for inflationary CDM 
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models. The core of the paper is Section 3, where we combine the COBE normahzation with 
the Lya P{k) to identify acceptable regions of the CDM parameter space. We focus on four 
representative models: a low density {^m < 1) open model, a low density flat model with a 
cosmological constant, and Einstein-de Sitter (0^ = 1) models with pure CDM and with a mixture 
of CDM and hot dark matter. Because different parameters have nearly degenerate influences 
on the predicted Lya P{k), we are able to summarize our results in terms of simple equations 
that constrain combinations of these parameters. In Section 4, we consider other observational 
constraints that can break these degeneracies, such as the cluster mass function, the peculiar 
velocity power spectrum, the shape of the galaxy power spectrum, and the CMB anisotropy power 
spectrum. We review our conclusions in Section 5. 



2. Parameter Space for CDM Models 



In simple inflationary models, the power spectrum of density fluctuations in the linear 
regime can be well approximated as a power law, P{k) oc A;" (where n = 1 is the scale-invariant 
spectrum), multiplied by the square of a transfer function T{k) that depends on the relative energy 
densities of components with different equations of state. We will assume the standard radiation 
background (microwave background photons and three species of light neutrinos) and consider as 
other possible components cold dark matter, baryons, a "cosmological constant" vacuum energy, 
and neutrinos with a non-zero rest mass in the few eV range. Within this class of models, the 
shape of the power spectrum is therefore determined by the parameters n, f^cDM) ^b, ^A, ^u, and 
h (since px = ^xPc oc fi^/i^). In place of Qb and rJcDM) we use the parameters 

B = n^h'^, (1) 

which is constrained by light-element abundances through big bang nucleosynthesis ( [Walker et al 



19911) , and 

= f^CDM + ^b + ^u, (2) 

which fixes r^cDM once B, h, and are specified. For non-zero 0,^, we assume one dominant 
family of massive neutrinos. We do not consider arbitrary combinations of and f^A but instead 
restrict our attention to the two theoretically simplest possibilities, spatially flat models with 
JIa = 1 — and open models with f^A = 0. 

Once the cosmological parameters are specifled, normalizing to the results of the COBE-DMR 
experiment determines the amplitude of P{k). For inflation models with n < 1, the COBE 
normalization can also be affected by the presence of tensor fluctuations (gravity waves). We 
consider normalizations with no tensor contribution and normalizations with the quadrupole 
tensor-to-scalar ratio T/ S = 7(1 — n) predicted by simple power law inflation models (e.g., 
Davis et al. 199^ ), but we do not consider arbitrary tensor contributions. We compute the 



COBE-normalized, linear theory, matter power spectrum P{k) using the convenient and accurate 
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fitting formulas of Eisenstein & Hu (1999), with the normahzation procedures of Bunn & White 
(1997) for all flat cases and for the open case without a tensor contribution and Hu & White 
(1997) for the open case with a tensor contribution. 

There are plausible variants of this family of inflationary CDM models that we do not analyze 
in this paper, because we lack the tools to easily calculate their predictions and because they 
would make our parameter space intractably larger. Prominent among these variants are models 
with a time-varying scalar field, a.k.a. "quintessence" (e.g., Peebles &: Ratra 1988 ; Wang fc 



Steinhardt 199^ ), models in which the energy of the radiation background has been boosted above 



its standard value by a decaying particle species, a.k.a. "rCDM" (e.g., [Bond &: Efstathiou 19911) , 



and models in which inflation produces a power spectrum with broken scale invariance (e.g., Kates 
et al. 1995| ). Given the observational evidence for a negative pressure component from Type la 



supernovae (Riess et al. 1998; Perlmutter et al. 1999), the quintessence family might be especially 
interesting to explore in future work. 

In sum, the free parameters of our family of cosmological models are Qm, h, n, B, fi^y, ^}\, 
and T/S. We allow h, n, B, and iljy to assume a continuous range of values. For and T/S 
we consider only two discrete options: Q\ = or 1 — ^Im, and r/5 = or 7(1 — n). 



3. Cosmological Parameters and the Lya Forest P(k) 



To organize our discussion and guide our analysis, we focus on variations about four fiducial 
models, each motivated by a combination of theoretical and observational considerations. The 
fiducial models are a flat cold dark matter model with a non-zero cosmological constant (ACDM), 
an open cold dark matter model with no cosmological constant (OCDM), an = 1 cold dark 
matter model with a significantly "tilted" inflationary spectrum (TCDM), and an Q^n = 1 model 
with a mixture of cold and hot dark matter (CHDM). 

For all of the fiducial models we adopt B = 0.02, based on recent measurements of the 
deuterium abundance in high-redshift Lyman limit absorbers ( Buries &: Tytler 1997] , 1998). For 
the TCDM and CHDM models we adopt h = 0.5 in order to obtain a reasonable age for the 
universe given the assumption that flm = 1- For the ACDM and OCDM models we instead adopt 
h = 0.65, which is better in line with recent direct estimates of the Hubble constant (e.g.. Mould 
|et al. 200q ). For the ACDM model we take Qrn = 0.4, but for OCDM we adopt a rather high 
density, = 0.55, in anticipation of our results in Section 4, where we consider the cluster mass 
function as an additional observational constraint. For the CHDM model, we take = 0.2 and 
assume one dominant species of massive neutrino; for all other models 17^^ — 0. ^Vith B^ /i, ^rnj 
and i^i, fixed, we are left with one free parameter, the inflationary spectral index n, which we 
choose in order to fit the amplitude of the Lya P(k) while maintaining the COBE normalization. 
The required value of n is different for models with no tensor contribution to CMB anisotropies 
than for models with tensor fluctuations; we refer to the fiducial models with tensor fluctuations as 
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ACDM2, 0CDM2, and TCDM2. Because a value n > 1 is required for CHDM and the assumption 
that T/ S = 7(1 — n) therefore cannot be correct in this case, we do not consider a CHDM model 
with tensor fluctuations. Table Q lists the parameters of the fiducial models. For later reference, 
Table |l] also lists each model's value of cts, the rms linear theory mass fluctuation in spheres of 
radius 8/i^^Mpc at z = 0. 

Figure ^ compares the power spectra of our fiducial models to the Lya P{k), shown as the 
filled circles with error bars. Note that the overall normalization of the data points is uncertain; 
at the la level they can shift up or down coherently by the amount indicated by the error bar on 
the open circle (see CWPHK for details). The COBE normalization itself has a la uncertainty 
of approximately 20% in P{k), roughly half of the Lya P{k) normalization uncertainty. Panels 
(a) and (c) show the fiducial models with and without tensors, respectively, over a wide range of 
wavenumber. Panels (b) and (d) focus on the range of wavenumbers probed by the Lya P{k). 
Our first major result is already evident from Figure ||: all of the fiducial models reproduce the 
observed Lya P{k). Each model has a single adjustable parameter, the spectral index n, so their 
success in reproducing both the amplitude and slope of P{k) is an important confirmation of a 
generic prediction of the inflationary CDM scenario, a point we will return to shortly. 

Within the precision and dynamic range of the CWPHK measurement, the Lya P{k) can be 
adequately described by a power law. CWPHK find 



with 




^\k)^—^P{k)=l^\k,)[-\ , (3) 



kp = 0.008 (km s"^ )-\ (4) 
/^\kp) = 0.573l°:?i, (5) 
V = -2.25 ±0.18. (6) 

Here A^(A;) is the contribution to the density variance per unit interval of Infc, and kp is a "pivot" 
wavenumber near the middle of the range probed by the data. 



Table 1: Fiducial Models 



Model 




Ha 


h 


n 


B 






T/S 


cs 


ACDM 


0.4 


0.6 


0.65 


0.96 


0.02 


0.0 







0.91 


ACDM2 


0.4 


0.6 


0.65 


0.98 


0.02 


0.0 


7(1 


-n) 


0.89 


OCDM 


0.55 


0.0 


0.65 


0.88 


0.02 


0.0 







0.67 


OCDM2 


0.55 


0.0 


0.65 


0.92 


0.02 


0.0 


7(1 


-n) 


0.64 


TCDM 


1.0 


0.0 


0.50 


0.84 


0.02 


0.0 







0.77 


TCDM2 


1.0 


0.0 


0.50 


0.89 


0.02 


0.0 


7(1 


-n) 


0.73 


CHDM 


1.0 


0.0 


0.50 


1.10 


0.02 


0.2 







0.96 
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Fig. 1. — Power spectra of the fiducial models (smooth curves) compared to the Lya P{k) 
determined by CWPHK (filled circles). Upper panels show models with no tensor fluctuations 
and lower panels show models with tensor fluctuations; the right hand panels zoom in on the 
range of wavenumbers probed by the Lya P{k). The error bar on the open circle shows an overall 
normalization uncertainty in the Lya P{k); at the la level all of the points can shift up or down 
coherently by this amount. Model parameters are listed in Table |. 
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In each panel of Figures ^ and ^, the central star shows the best fit values of A^{kp) and v 
quoted above, and the two large concentric circles show the la (68%) and 2a (95%) confidence 
contours on the parameter values. The calculation of these confidence contours is described in 
detail in Section 5 of CWPHK. Briefly, the likelihood distribution for the slope, u, is derived by 
fitting the power law form (eq. ^) to the P{k) data points, using their covariance matrix. The 
likelihood distribution for the amplitude, A'^{kp), is obtained by convolving the distributions 
calculated from two separate sources of uncertainty involved in the P{k) normalization. The 
joint confidence contours on the two parameters are obtained by multiplying together the two 
independent likelihood distributions. The la and 2a contours correspond to changes in — 21ogg£ 
from its best fit value of 2.30 and 6.17, respectively, where C is the likelihood. 

The open circular point near the middle of each panel of these figures shows the fiducial 
model's prediction of A'^{kp) and u. ACDM, OCDM, and TCDM models without tensors appear 
in the left column of Figure ^, the corresponding models with tensors appear in the right column of 
Figure ^, and the CHDM model appears in Figure |3|. As expected from Figure |l|, the fiducial model 
predictions lie well within the 68% confidence contour in all cases. The 20% COBE normalization 
uncertainty adds a log(1.2) ~ 0.08 error bar to the predicted value of log A^(A;p), which we have 
not included in the plots. Because this uncertainty is small (once added in quadrature) compared 
to the Lya P{k) uncertainty itself, we have ignored it in the analysis of this paper. With a higher 
precision Lya forest measurement, it would be important to include the COBE normalization 
uncertainty as an additional source of statistical error. 

Changing any of the parameter values in any of the models shifts the predicted A'^{kp) and u, 
and the remaining points in Figures |2| and ^ show the effects of such parameter changes. Taking 
the ACDM model of Figure ^ as an example, the two filled circles show the effect of increasing 

by 0.1 and 0.2 (to J7m = 0.5 and = 0.6), while maintaining the condition + ^^a = 1 and 
keeping all other parameters fixed at the fiducial values listed in Table ||. The two open circles 
show the effect of decreasing by 0.1 and 0.2. With = 0.2 and other parameters unchanged 
(leftmost open circle), the predicted amplitude A^(fcp) falls below the 95% confidence lower limit 
of CWPHK. In similar fashion, filled (open) pentagons show the effect of increasing (decreasing) 
n by 0.05, filled (open) squares show the effect of increasing (decreasing) h by 0.05, and filled 
(open) triangles show the effect of increasing (decreasing) by 0.01, in all cases keeping the other 
parameters fixed at their fiducial values. The format of the other panels of Figure ^ is identical, 
except that we do not show 0^ changes for TCDM. For ACDM2 and 0CDM2, we do not allow 
n > 1. In Figure |3|, filled (open) hexagons show the effect of increasing (decreasing) 0.,^ by 0.1 
while keeping Qrn = 1- Open circles show the effect of decreasing fi^ by 0.1 while adding Q\ to 
maintain fiat space; results are virtually indistinguishable if 0,\ is zero and the universe becomes 
(slightly) open. We do not consider changes that make ^Im > 1- 

Parameter changes have similar effects in all of the models, and these effects can be easily 
understood by considering the physics that determines the shape and normalization of the matter 
power spectrum. The CDM transfer function has a single fundamental scale cteq determined 




-0.8 -0.6 -0.4 -0.2 0.2 -0.8 -0.6 -0.4 -0.2 0.2 
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-2.5 



-0.8 -0.6 -0.4 -0.2 0.2 -0.8 -0.6 -0.4 -0.2 0.2 
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-2.5 
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log A=(kJ log A=(k^) 

Fig. 2. — Constraints on the parameters of CDM models from COBE and the Lya P{k) 
measurement. In each panel, the central star shows CWPHK's best-fit values of A^{kp) and 
and closed contours show the 68% and 95% confidence regions. Each panel corresponds to a 
different one of the fiducial models, with the central open circle marking the model prediction 
for the parameters listed in Table |l[ Other filled (open) points show the effects of increasing 
(decreasing) these parameters by fixed amounts while keeping all other parameters fixed. Circles 
show changes Ailm = 0.1, pentagons An = 0.05, squares A/i = 0.05, and triangles Afi^ = 0.01; 
f^m changes are not considered for TCDM, and n > 1 is not considered for models with tensor 
fluctuations. The error cross shows the Lya P{k) measurement of McDonald et al. (2000). 
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Fig. 3.- 



Like Figure |2|, but for the CHDM model. Hexagons show changes of by Afi^ = 0.1. 
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by the size of the horizon at the time of matter-radiation equahty; this scale is roughly the 
wavelength at which the power spectrum turns over. Increasing h, and hence the matter density 
Pm oc Q^h'^i moves matter-radiation equality to higher redshift and lower igqi shifting the model 
power spectrum towards smaller scales (to the right in Figure |^). This horizontal shift, combined 
with an upward vertical shift to maintain the COBE normalization on large scales, increases the 
amplitude of P{k) on Lya forest scales and translates a shallower (higher u) part of the spectrum 
to kp. Increasing also lowers teq and therefore has a similar effect. Open models are more 
sensitive than flat models to changes in because the integrated Sachs- Wolfe effect makes a 
greater contribution to large scale CMB anisotropies ( ^achs fc Wolfe 1967 ; Hu, Sugiyama, fc Silk 



1997). Increasing reduces the integrated Sachs- Wolfe effect and hence increases the matter 



fluctuation amplitude implied by COBE, shifting the power spectrum vertically upward. The 
value of A'^(kp) is sensitive to the spectral index n because of the very long lever arm between 
the COBE normalization scale and the scale of the Lya forest measurement. A small decrease 
in n produces an equally small decrease in v but a large decrease in A'^{kp). The fluctuation 
amplitude is even more sensitive to n in tensor models because, with T/S = 7(1 — n), decreasing n 
also increases the contribution of gravity waves to the observed COBE anisotropies and therefore 
reduces the implied amplitude of the (scalar) matter fluctuations. Since fluctuations in the baryon 
component can only grow after the baryons decouple from the photons, increasing B depresses 
and steepens P{k) on small scales and therefore reduces A'^{kp) and u. However, for our adopted 
parameters the baryons always contribute a small fraction of the overall mass density, so the 
influence of Of, changes is small. Increasing Qi, in the CHDM model has a much greater effect in 
the same direction, since the suppression of small scale power by neutrino free streaming is much 
greater than the suppression by baryon-photon coupling. 

Figures |2| and |^ re-emphasize the point made earlier in our discussion of Figure the 
agreement between the predicted and measured slope of the Lya P{k) confirms a general prediction 
of the inflationary CDM scenario. Although the four fiducial models correspond to quite different 
versions of this scenario, all of them reproduce the measured value of = —2.25 to well within 
its la uncertainty once the value of n is chosen to match the measured A'^{kp). However, if 
the measured value of v had been substantially different, e.g. implying i/ > —2 or < —2.5, 
then none of these models could have reproduced the measured while remaining consistent 
with the measured A'^{kp), even allowing chang es in 72, ^^772, ^i/, or r^f,. A different value of z/ 
would therefore have been a challenge to the inflationary CDM scenario itself rather than to any 
specific version of it. Note also that any of the models would match the observed v within its la 
uncertainty even if we had assumed a scale- invariant, n = 1 inflationary spectrum; it is the A'^{kp) 
measurement that requires the departures from n = 1. Because of the long lever arm from COBE 
to the Lya P{k), parameter changes that have a modest effect on u have a large effect on A'^{kp). 

Figure ^ also shows that changes of the different model parameters have nearly degenerate 
effects on the predicted values of A'^{kp) and u. For example, in the ACDM model, increasing 0^ 
by 0.1 would increase the predicted slope and amplitude, but decreasing h by 0.05 would almost 
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exactly cancel this change. This near degeneracy allows us to summarize the constraints imposed 
by COBE and the Lya P{k) with simple formulas of the form 

flmh'^n^B^ = k±e, (7) 

where k is the value obtained for the best-fit parameter values in Table |l] and the uncertainty e 
defines the variation that is allowed before the model leaves the 68% confidence contour. Table Q 
lists the values of a, P, 7, k, and e for all of the fiducial models. Although we do not show 
changes for the TCDM models in Figure |2|, we vary it below 1.0 (adding to keep the universe 
flat) in order to derive the a, (5, and 7 indices, so that in all models their values reflect the 
importance of a change in /i, n, or B relative to a change in Jim- 
Equation together with Table |2|, is our second principal result, defining the quantitative 
constraints placed on the parameters of inflationary CDM models by the combination of COBE 
and the Lya forest P{k). The values of the a, /?, and 7 indices reflect the sensitivity of the 
predicted power spectrum amplitude A^(/cp) to the model parameters, quantifying the impressions 
from Figure |2|. Again taking ACDM as an example, we see that small variations in h and n have 
much greater effect than small variations in Vim, and that the suppression of small scale power 
from increases in B is always a modest effect. Models with tensors are much more sensitive to n 
than models without tensors because of the influence of gravity waves on the P{k) normalization, 
as discussed above. Although the index values are derived in all cases by considering small 
variations about the corresponding fiducial model, the constraint formula ^ remains accurate 
even for fairly large changes in the cosmological parameters. For example, plugging the TCDM 
values of Q.m, h, n, B into equation (0) with the ACDM values of a, f3, and 7 yields k = 0.47, 
compared to the value k = 0.44 listed for ACDM in Table ^. 

Figure ^ shows that the effects of parameter changes are less degenerate in the CHDM 
model. This difference in behavior is not surprising, since neutrino free streaming changes P{k) by 
depressing it at small scales rather than simply shifting or tilting it. The slope v is therefore much 
more sensitive to changes in Qi, than to changes in other parameters. Conversely, the influence of 
h on 1/ through shifting tgq is nearly cancelled by the effect of h on the implied neutrino mass and 



Table 2: Constraint Parameters (see equation ^ 



Model 


a 


/3 


7 


S 


k 


€ 


ACDM 


1.88 


2.68 


-0.26 




0.44 


0.12 


ACDM2 


1.84 


4.48 


-0.25 




0.43 


0.08 


OCDM 


1.55 


2.57 


-0.23 




0.50 


0.10 


OCDM2 


1.80 


3.45 


-0.18 




0.38 


0.08 


TCDM 


2.33 


4.08 


-0.40 




0.46 


0.14 


TCDM2 


1.82 


4.60 


-0.15 




0.30 


0.08 


CHDM 


0.93 


1.74 


-0.13 


-0.37 


1.87 


0.26 
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free streaming length. We still analyze this case as above, adding a factor of , to obtain 

n^h^n^B^n^^ = k±e, (8) 

with parameters also listed in Table |^. However, this equation cannot describe the results of 
Figure |3| as accurately as equation ^ describes the results of Figure 

Recently McDonald et al. (2000) measured the Lya forest flux power spectrum in a sample 
of eight Keck HIRES spectra and used it to infer the amplitude and shape of the mass power 
spectrum. Their mean absorption redshift is z ~ 3 rather than z = 2.5, and their data best 
constrain the P{k) amplitude at A; = 0.04 (km s^^ rather than 0.008 (km s^^ )^^. However, 
assuming gravitational instability and a CDM power spectrum shape, they extrapolate from their 
result to derive values of v and that can be directly compared to CWPHK's measurement 
at z = 2.5, kp = 0.008 (km s"^ )-\ obtaining u = -2.24 ± 0.10 and A'^{kp) = 0.32 ± 0.07. 
Despite the entirely independent data sets and very different modelling procedures, the CWPHK 
and McDonald et al. (2000) measurements agree almost perfectly in slope and are consistent 
in amplitude at the ~ la level. We plot the McDonald et al. (2000) measurement as error 
crosses in Figures || and ^. McDonald et al. (2000) note that the small error bar on should 
be considered preliminary, since they have not fully investigated the sensitivity of their power 
spectrum normalization procedure to their modelling assumptions. 

Clearly none of our qualitative conclusions about inflationary CDM models would change 
if we were to adopt the McDonald et al. (2000) P{k) determination instead of the CWPHK 
determination. Conveniently, the McDonald et al. (2000) point lies almost exactly on our —la 
error contour, so to a good approximation one can obtain the parameter constraints (0) and (|8|) 
implied by the McDonald et al. (2000) measurement by simply replacing the values of k in Table ^ 
by A; — e. 



4. Combining with other constraints 



We have shown that the combination of COBE and the Lya P{k) yields constraints 
on degenerate combinations of cosmological parameters. To break these degeneracies, we now 
consider observational constraints from other studies of large scale structure and CMB anisotropies. 
Analyses of cosmological parameter constraints from multiple observations have been carried out 



by numerous groups (recent examples include Bahcall et al. 1999; Bridle et al. 1999; Steigman 



Hata, fc Felten 1999| ; [NovosyadlyJ et al. 2000| ; [Wang, Tegmark, fc Zaldarriaga 200T|) . Our new 
contribution is to include the Lya ^(A;) as one of the observational constraints (also considered 
by [Novosyadlyj et al. 200C and Wang, Tegmark, fc Zaldarriaga 2001)) . We focus our attention on 
several other constraints that can be cast into a form that complements our results from Section 
3: the mass function of galaxy clusters, the mass power spectrum inferred from galaxy peculiar 
velocities, the shape parameter of the galaxy power spectrum, and a constraint on n from CMB 
anisotropy data. Our discussion in this Section will be more qualitative than our discussion in 
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Section 3, in part because the uncertainties in these constraints are largely systematic, so that a 
straightforward statistical combination could be misleading. 

In each panel of Figures ^ and ^, the heavy solid line shows the locus of {Qm, n) values that 
yield a simultaneous match to COBE and the CWPHK measurement of the Lya P{k). These 
lines are very close to those implied by equation ^ and Table but since those results are, 
strictly speaking, expansions about our fiducial model parameters, we compute the best-fit value 
of n exactly for each rather than using equation (0). The ztlo" constraints are shown as 
the lighter solid lines; these are close to the curves implied by equation ^ and Table ^ with k 
replaced by A; it e. Because the Lya P{k) constraint is not very sensitive to we keep B fixed 
at our fiducial value of 0.02 in all cases. We show results for h = 0.65, h = 0.45, and h = 0.85 in 
the upper, middle, and lower panels of each figure, with fiat and open models in the left and right 
hand columns, respectively. Figure ^ shows models without tensors and Figure |^ models with 
tensors. For models with tensors, we restrict the parameter space to n < 1, since our assumption 
that T/ S = 7(1 — n) only makes sense in this regime. The TCDM models can be considered as 
the limit of either the fiat or open models at = 1- Note that the McDonald et al. (2000) 
estimate of the Lya P{k) corresponds very closely to our —la constraint, so to adopt McDonald 
et al. (2000) instead of CWPHK one simply follows the lower solid line instead of the middle solid 
line as the constraint. 

For Gaussian initial conditions, the space density of clusters as a function of virial mass 
constrains a combination of Q.m and the mass fiuctuation amplitude, since clusters of a given mass 
can be formed by the collapse of large volumes in a low density universe or smaller volumes in a 
higher density universe. This constraint can be summarized quite accurately in a formula relating 
Q^rn to the rms mass fluctuation ( [White, Efstathiou, fc Frenk 1993a| ). We use the specific version 
of this formula obtained by Eke, Cole, & Frenk (1996, hereafter ECF) using N-body simulations 
and the Press-Schechter (1974) approximation: 

as = (0.52 ± 0.04)O^-°-^6+°-i°^™ = 

as = (0.52 ± 0.04)O^-0-52+0-i3f^m J^A = 1 - . ^' 

For each value of 17^, we find the value of ug required by the cluster mass function from 
equation (^). Given h and B = 0.02, we then find the value of n required to produce this value of 
(Tg by numerically integrating the CDM power spectrum. This constraint in the — n plane is 
shown by the dotted line in each panel of Figures ^ and |5|, with an error bar that indicates the 8% 
uncertainty quoted in equation @) from ECF. 

For a given value of f^mi the matter power spectrum can also be estimated from the statistics 
of galaxy peculiar motions. Freudling et al. (1998) apply a maximum likelihood technique to the 
SFI peculiar velocity catalog to constrain COBE-normalized, infiationary CDM models for the 
matter power spectrum, obtaining the constraint 
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Fig. 4. — Constraints in the — n plane from a variety of cosmological tests, for models with 
no tensor fluctuations. Upper, middle, and lower panels show models with h = 0.65, h = 0.45, 
and h = 0.85, respectively, and in all cases we keep B = 0.02. Flat models appear in the left 
hand column, open models with = in the right hand column. In each panel, the heavy solid 
line shows the 0^ — n locus determined by the combination of COBE and the Lya P{k)^ and 
the light solid lines show the itlo" range of this locus. Dotted lines show the constraint (P) from 
the cluster mass function, short-dashed lines the constraint ( |lO|) from the peculiar velocity power 
spectrum, long-dashed lines the constraint (lTl| ) from the shape of the galaxy power spectrum, and 
horizontal dot-dashed lines the constraint on n from CMB anisotropy measurements. Error bars 
show representative la statistical uncertainties in these constraints. A model is consistent with 
multiple constraints if it lies in the region of the Vim — n plane where these constraints overlap 
within their uncertainties. 
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Fig. 5. — Like Figure ^ but for models with tensor fluctuations. Note that there is no version of 
the velocity power spectrum constraint for 0CDM2. The CMB constraint on n does not impose the 
condition T/S = 7(1 — n), but all of the other constraints do assume this condition and therefore 
cannot be applied for n > 1. 
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where /i, ly, k and e are dependent on the cosmology and /igo = h/O.Q. In a flat, Q\ = 1 — 
model with no tensor component, (/x, u, k, e)=(1.3, 2.0, 0.58, 0.08), while if a tensor component is 
allowed they become (1.3, 3.9, 0.58, 0.08). For an open, = model without a tensor component 
they are (0.9, 1.4, 0.68, 0.07). Freudling et al. (1998) do not consider open, JIa = cases with a 
tensor component. For specified h, equation (^) yields a constraint in the Qm — n plane, shown 
by the short-dashed line in the panels of Figures ^ and |5|. The associated la error bars are based 
on the statistical uncertainties e quoted by Freudling et al. (1998) and listed above. For brevity, 
we will refer to these curves as the velocity power spectrum constraint, though they represent the 
constraints on the density power spectrum implied by peculiar velocites. 

We do not want to use the amplitude of the galaxy power spectrum as one of our constraints 
because it can be strongly affected by biased galaxy formation. However, a variety of analytic 



and numerical arguments (e.g., Coles 1993; 


Fry & Gaztahaga 1993; 


Mann, Peacock, & Heaven^ 


1998 


; Scherrer & Weinberg 1998 


; Narayanan, Berlind, & Weinberg 2000 


) suggest that biased 



galaxy formation should not alter the shape of the galaxy power spectrum on scales in the linear 
regime, and on these scales the shape is directly related to the parameters of the inflationary 
CDM cosmology. We adopt the specific constraint found by Peacock & Dodds (1994) from their 
combined analysis of a number of galaxy clustering data sets: 



-0 J 1 + 



2h 



0.32 



1 



n 



1 



0.255 ±0.017. 



(11) 



For n = 1, Fcfj = F, where F is the shape parameter in the conventional parameterization of 
the inflationary CDM power spectrum ( Pardeen et al. 1986 ; the influence of Qb is discussed 
by ^ugiyama 1995 ). While the effects of F and n on the power spectrum shape are different. 



equation (11) combines them in a way that approximates their nearly degenerate influence over 
the range of scales currently probed by large scale clustering measurements. For specified h and B, 
equation ([lT| ) becomes a constraint in the Qm — n plane. We plot this constraint as the long-dashed 
line and associated error bar in the panels of Figures ^ and ||. We should note, however, that 
the Peacock &: Dodds (1994) error bar may be overoptimistic, since independent estimates of 
Feff often fall outside this range. Eisenstein & Zaldarriaga (2001) have recently re-examined the 
spatial power spectrum inferred from the APM survey and conclude that the 68% confidence 
interval of F (for n = 1) is 0.19 — 0.37, much larger than the range implied by equation ([TT|), and 
Efstathiou & Moody (2000) favor a lower central value (F 0.12, with 2a range 0.05 < F < 0.38). 
As older estimates of the galaxy power spectrum are supplanted by results from the 2dF and 
Sloan galaxy redshift surveys, the F parameterization itself may become an insufficiently accurate 



representation of the theoretical predictions ( Percival et al. 2001 ) 



A detailed consideration of constraints from smaller scale CMB anisotropy measurements is 
beyond the scope of this paper, but we do want to draw on limits that smaller scale measurements 
place on the inflationary index n. For the no-tensor models, we adopt the "weak prior" constraint 



n 



0.96to:o9 of Netterfield et al. (2001), based on data from the BOOMERANG experiment. 



which we represent by the horizontal dot-dash line and la error bar in Figure 0. Since Netterfield 
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et al. (2001) do not consider models with tensor fluctuations, we take the corresponding constraint 
for the tensor models in Figure ^ from Wang et al. (2001). Their model space is less restrictive 
than ours because they do not impose the power-law inflation relation T/ S = 7(1 — n), and using 
CMB data alone they find only a very weak constraint on n. We therefore adopt their constraint 
from the combination of CMB and large scale structure data, n = 0.911q'q5, where we have 
reduced the 95% confidence range quoted in their table 2 by a factor of two to get a representative 
Ifj uncertainty. 

In Figures |^ and ^, the cluster mass function, velocity power spectrum, and shape parameter 
constraints tend to be roughly parallel to each other, with the shape parameter following a 
somewhat different track when tensor fluctuations are important. The shape parameter constraint 
is usually compatible with the cluster mass function constraint, at least if one allows for the 
possibility that the error bar in equation (|l^) is somewhat too small. However, the velocity power 
spectrum always implies a higher fluctuation amplitude than the cluster mass function, and the 
two constraints are not consistent within their stated la uncertainties for any combination of fi^, 
n, and h. A recent analysis by Silberman et al. (2001) shows that the discrepancy is probably a 
result of non-linear effects on the velocity power spectrum, and that correcting for these yields 
results closer to the cluster constraint. We therefore regard the cluster constraint as more reliable, 
and we retain the velocity power spectrum curve mainly as a reminder of other data that can be 
brought to bear on these questions. 

The Lya P{k) curve cuts across the other three constraints, requiring greater change in 0.^ 
for a given change in n. The CMB anisotropy limit on n cuts across all of the other constraints. 
The COBE-DMR measurement is represented implicitly in Figures ^ and |5| through its role in the 
Lya P{k) constraint, the velocity power spectrum constraint, and the CMB anisotropy constraint 
on n. The size of the la error bars in these flgures, and the probability that at least some of them 
are underestimated, prevents us from drawing sweeping conclusions. However, Figures I and I 
do have a number of suggestive implications if we look for models that lie within the overlapping 
Ifj uncertainties of the various constraints. Since it is not possible to satisfy the cluster mass 
function and velocity power spectrum constraints simultaneously within the class of models that 
we consider, the implications depend strongly on which of these constraints we take to be more 
reliable. The shape parameter implications are usually intermediate, but significantly closer to 
those of the cluster mass function. 

The combination of the velocity power spectrum and Lya P{k) constraints implies a high 
density universe, with ^ 1 preferred and Vl^ ^ 0.6 separating the two constraints by more 
than their la error bars. The Lya P{h) constraint rules out the high values of n that could 
otherwise allow low in equation ([lO|). For h > 0.65, intersection of the velocity power spectrum 
and Lya P{k) constraints occurs at n ^ 0.8, incompatible with the CMB anisotropy constraint. 
However, an ~ 1 universe would require a low value of h in any case because of the age 
constraint for globular cluster stars, and this would push the intersection to higher n. As noted 
earlier, the velocity power spectrum constraint shown here is probably biased towards high by 
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the non-linear effects described by Silberman et al. (2001). 

If we instead adopt the cluster mass function constraint, then consistency with the Lya P{k) 
and COBE requires < 1- For h = 0.65, the constraints intersect at ^Im ~ 0.4 — 0.5 in flat 
models and ~ 0.5 — 0.6 in open models; increasing h slightly decreases the preferred O-m and 
vice versa. This conclusion — that the combination of COBE, the Lya -P(fc), and the cluster mass 
function implies a low density universe — is the most important and robust result to emerge from 
this multi-constraint analysis. 

At one level, our conclusions about the matter density come as no suprise, since we have 
already argued, in Weinberg et al. (1999), that consistency between the cluster mass function and 
the Lya -P(fc) implies 0,m in this range independent of the COBE normalization. However, the 
nature of the argument is subtly different in this case. In Weinberg et al. (1999), we considered 
matter power spectra of the CDM form parameterized by T (with n = 1), and by combining the 
Lya -P(fc) measurement with the cluster constraint (^), we found = 0.34 -|- 1.3(r — 0.2) for 
flat models and = 0.46 -|- 1.3(r — 0.2) for open models, with la uncertainties of about 0.1. 
However, the Lya P{k) alone could not rule out the solution of high and high F, so Weinberg 
et al.'s (1999) conclusion that < 1 rested crucially on the empirical evidence for F ~ 0.2 from 
the shape of the galaxy power spectrum. Within the class of CDM models considered here, the 
combination of COBE and the Lya P{k) determines n, and hence the effective value of F (eq. |Tll) , 
once Qm, h, and B are specified. Simultaneous consistency between COBE, the Lya P{k), and the 
cluster mass function requires low independent of the galaxy power spectrum shape, thereby 
strengthening the overall argument for a low density universe, and, by the by, for a matter power 
spectrum with low T^s- The lower limit on from this combination of constraints varies with 
the choice of other parameters, but it never reaches as low as flm = 0.2 unless h > 0.85. 

For all of the models shown in Figures Q and ||, the Lya P{k) and cluster mass function 
constraints intersect at values of n consistent with the CMB anisotropy constraints, provided one 
takes the 1-a error ranges into account. A factor of two improvement in the precision of the Lya 
P{k) measurement could greatly restrict the range of models compatible with all three constraints, 
especially if the Lya P{k) amplitude is somewhat lower, as McDonald et al. (2000) find. 

There are, of course, numerous other constraints on cosmological parameters, and we will 
briefly consider three of them: the cluster baryon fraction, the location of the first acoustic peak in 
the CMB power spectrum, and the evidence for accelerating expansion from Type la supernovae. 
(Our focus on h = 0.65 as a fiducial case already reflects our assessment of the most convincing 
direct estimates of Hq.) If one assumes that baryons are not overrepresented relative to their 
universal value within the virial radii of rich clusters, then the combination of the measured gas 
mass fractions with big bang nucleosynthesis limits on yields an upper limit on $7^ ( [White et 
al. 1993b). Applying this argument, Evrard (1997) concludes that 



nn^nr^h-'^/^ < 23.1 =^ n^< 0.57 ( (J^) ' , (12) 



0.02/ V0.65 



- 19 - 



at the 95% confidence level. From Figures ^ and ^ we see that models matching COBE, the Lyct 
P{k), and the cluster mass function are always consistent with this limit — easily in the case of 
flat models, sometimes marginally in the case of open models. Models that match the velocity 
power spectrum instead of the cluster mass function are usually incompatible with this limit, 
though sometimes only marginally so. 

The location of the first acoustic peak in the CMB anisotropy spectrum is a strong diagnostic 
for space curvature (e.g., Doroshkevich, Zeldovich, &: Sunyaev 1978| ; [Wilson fc Silk 1981 ; Sugiyama 



fc Gouda 1992| ; Kamionkowski, Spergel, fc Sugiyama 1994 ; Hu et al. 1997), and recent anisotropy 



measurements on degree scales favor a geometry that is close to flat (e.g., [Miller et al. 1999 



Melchiorri et al. 2000 ; Netterfield et al. 2001 ; Pryke et al. 2001 ). Clearly our flat universe models 



are compatible with these results, as are the open universe models that match the Lya P{k) and 
the velocity power spectrum (all of which have close to one). The open models that match 
Lya P{k) and the cluster mass function are generally ruled out by the most recent, high precision 
limits on space curvature. The Type la supernova measurements of the cosmic expansion history 
(Riess et al. 1998; Perlmutter et al. 1999) add a great deal of discriminatory power, since they 



constrain a parameter combination that is roughly — instead of + ^^a; Perlmutter 
et al. (1999) quote Qm — O.TSf^A ~ —0.25 it 0.125. All of the open models miss this constraint 
by many a, and the flat models matching the Lya P{k) and the velocity power spectrum fail 
because the values of Qm are too high. The combination of COBE, the Lya P{k), and the cluster 
mass function, on the other hand, is compatible with the supernova results for flat models with a 
cosmological constant, though it favors somewhat higher values of 

We have not carried out a similar multi-constraint analysis for the CHDM model because the 
formulas ( [Til ) and (|^) for the shape parameter and velocity power spectrum constraints do not 
apply to it and the formula (^) for the cluster mass function constraint may be less accurate for 
non-zero il.^- However, our fiducial CHDM model, with = 0.2, has as = 0.96, with n = 1.10. 
For Qi, = 0.3 we obtain as = 1.15 (n = 1.23), for Qi, = 0.1 we obtain as = 0.81 (n = 0.96), and for 
the TCDM model, which represents the limiting case of ^li, = 0, we obtain as = 0.77 (n = 0.84). 
All of these models are likely to violate the cluster mass function constraint, which according to 
equation (|) imphes as = 0.52 ± 0.04 for = 1- We conclude that COBE-normahzed CHDM 
models with = 1, /i ~ 0.5 cannot simultaneously match the Lya P{k) and the cluster mass 
function. The Lya P{k) strengthens the case against this class of CHDM models by ruling out the 
low values of n that would otherwise allow them to match cluster masses ( [Ma 1996[) . Of course 



CHDM models with $7^ < 1 can satisfy the observational constraints for appropriate parameter 
choices, and the general problem of using CMB measurements and the Lya P{k) to measure Qi, is 
discussed by Croft, Hu, & Dave (1999). However, the possible presence of a neutrino component 
does not alter our conclusion that COBE, the Lya P{k), and the cluster mass function together 
require a low density universe. 

All in all, the CWPHK and McDonald et al. (2000) measurements of the Lya P{k) provide 
additional support for the current "consensus" model of structure formation, ACDM with ~ 0.4 
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Fig. 6. — CMB anisotropy power spectra for five of the fiducial models that match COBE-DMR 
and the Lya P{k). Ci is the mean-squared amplitude of spherical harmonics of order /. The 
models shown are ACDM (solid), ACDM2 (dotted), OCDM (long-dashed), TCDM (short-dashed), 
and CHDM (dot-dashed). Parameters of the models are listed in Table |l[ Data points with Icr 
errors are taken from Wang et al. (2001). We shift from logarithmic spacing to linear spacing at 
/ = 100 in order to show both large and small angular scales clearly. 
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and h « 0.65. Moderate improvements in the statistical precision of the constraints considered 
here could strengthen this support, or they could open fissures of disagreement. Improvements in 
the near future could also allow some interesting new tests, such as discriminating between models 
with no tensor fluctuations and models with the T/S = 7(1 — n) contribution predicted by power 
law inflation. 

A detailed consideration of the constraints from CMB anisotropy measurements is a major 
undertaking in itself, well beyond the scope of this paper. However, to illustrate the interplay 
between our results and recent CMB experiments, we plot in Figure ^ the predicted CMB 
power spectra of five of our fiducial models: ACDM, ACDM2, OCDM, TCDM, and CHDM. We 
computed these power spectra using CMBFAST (ISeljak fc Zaldarriaga 1996t [Zaldarriaga, Seljak,| 



fc Bertschinger 199^ ), with the cosmological parameter values listed in Table ||. The CHDM 



model stands out from the rest because matching the Lya P{k) requires a high value of n, which 
boosts the anisotropy on small scales. The OCDM model also stands out, albeit less dramatically, 
because the open space geometry shifts the acoustic peaks to smaller angles. The TCDM model 
lies below the ACDM models because of its larger tilt, which suppresses small scale fluctuations. 
Figure |6| shows data points taken from the joint analysis of numerous CMB data sets by Wang et 
al. (2001; see their Table 1). The two ACDM models fit these data points remarkably well, given 
that the choice of their parameters was not based on small scale CMB data at all. Because the 
combination of COBE and the Lya P{k) implies n close to one for both of these models, their 
predictions are not very different, and the current CMB data do not distinguish between them. 
However, the TCDM, OCDM, and CHDM models are clearly ruled out, and while we have not 
attempted to adjust their parameters within the constraints allowed by equations (Q) and (^), it 
appears unlikely that any such adjustment would allow these models to fit the current CMB data. 



5. Conclusions 



The slope of the mass power spectrum inferred by CWPHK from the Lya forest, 
v = —2.25 lb 0.18 at kp = 0.008 (km s^^ at z = 2.5, confirms one of the basic predictions of the 
inflationary CDM scenario: an approximately scale-invariant spectrum of primeval inflationary 
fluctuations (n 1) modulated by a transfer function that bends the power spectrum towards 
P{k) oc k^~^ on small scales. If the measured slope of the power spectrum had implied u > —2 
or v < —2.5, we would have been unable to reproduce the Lya P{k) with any of the models 
considered here, even allowing wide variations in the cosmological parameters. 

Because the amplitude of the COBE-normalized power spectrum on small scales is very 
sensitive to n, we are able to match the CWPHK measurement of A^(fcp) in most of the major 
variants of the CDM scenario (ACDM, OCDM, TCDM, CHDM) by treating n as a free parameter. 
Within each of these variants, we obtain constraints on the model parameters of the form 
Qrnh'^n^ B'^ = k ziz e (eq. 0) or Q.rnh"n^ B'^'^y^ = k ziz e (eq. ^), with the parameter values listed in 
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Table |2|. These constraints, together with the confirmation of the predicted slope, are the main 
results to emerge from combining the Lya P{k) measurement with the COBE-DMR result. 

As shown in Figures ^ and ^, the parameter combination constrained by COBE and the 
Lya P{k) is different from the combinations constrained by other measurements of large scale 
structure and CMB anisotropy, so joint consideration of these constraints can break some of the 
degeneracies among the fundamental parameters. If we combine the Lya P{k) constraint with the 
constraint on 0^ and dg inferred from the cluster mass function (White et al. 1993a; ECF), then 
we favor a low density universe, with 0,rn ~ 0.3 — 0.5 in flat models and 0,m ~ 0.5 — 0.6 in open 
models. This combination is also consistent with CMB anisotropy constraints on n. The open 
models are inconsistent with the angular location of the first acoustic peak in the CMB power 



spectrum (Netterfield et al. 2001; Pryke et al. 2001), and they are strongly inconsistent with Type 



la supernova results, which imply - 0.75Oa « -0.25 ± 0.125 (p less et al. 1998| , |Perlmuttei] 



et al. 199^ ). The flat models are consistent with both constraints. On the whole, the CWPHK 



measurement of the Lya P{k) supports the consensus in favor of ACDM with ~ 0.4, h ~ 0.65. 
The contribution of the Lya P{k) to this consensus comes both from the slope, which confirms the 
generic inflationary CDM prediction, and from the amplitude, which has a different dependence 
on cosmological parameters than any of the other constraints considered here. 

There are bright prospects for improvements of this approach in the near future. McDonald 
et al. (2000) have inferred the mass power spectrum from an independent Lya forest data set 
using a different analysis method, obtaining a nearly identical slope and an amplitude lower 
by ~ la. We have recently analyzed a much larger data set of high and moderate resolution 
spectra, using a variant of the Croft et al. (1998, 1999) method, and the improved data yield 
much higher statistical precision and better tests for systematic effects. Constraints from this 
new measurement of P{k), using the method developed here, are presented in §7 of Croft et al. 
(2001). Recent measurements of CMB anisotropy have greatly improved the level of precision on 
small angular scales, and results from the MAP satellite should provide another major advance in 
the near future. These measurements yield tighter cosmological parameter constraints on their 
own, but they become substantially more powerful when combined with data that constrain the 
shape and amplitude of the matter power spectrum. It is evident from Figures ^ and ^ that 
simply reducing the error bars on n and the Lya P{k) by a factor of two would already produce 
interesting new restrictions on the allowable range of models. These restrictions can become 
very powerful if ongoing studies of cluster masses using galaxy dynamics. X-ray properties, the 
Sunyaev-Zel'dovich effect, and gravitational lensing confirm the robustness of the cluster mass 
function constraint. In the slightly longer term, the 2dF and Sloan redshift surveys should produce 
measurements of the shape of the galaxy power spectrum that shrink the current statistical and 
systematic uncertainties, so that demanding consistency between the inferred value of Fefr and 
other constraints becomes a useful additional test. At the very least, these developments should 
lead to a powerful test of the infiationary CDM picture and high-precision determinations of its 
parameters. If we are lucky, improved measurements will reveal deficiencies of the simplest ACDM 



-23- 



models that are hidden within the current uncertainties, and resolving these discrepancies will 
lead us to a better understanding of the cosmic energy contents and the origin of primordial 
fluctuations in the hot early universe. 
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